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* By N. Fuchs 


SUMMARY 

The evaporation velocity of liquid droplets under various conditions 
is theoretically calculated and a number of factors are investigated 
which are neglected in carrying out the fundamental equation of Maxwell'..,. 
It is shown that the effect of these factors at the small drop sizes and 
the small weight concentrations ordinarily occurring in fog can be calcu- 
lated by simple corrections. The evaporation process can be regarded as 
quasi-stationary in most cases. 

The question at hand, and also the equivalent question of the veloc- 
ity of growth of droplets in a supersaturated atmosphere, is highly sig- 
nificant in meteorology and for certain industrial purposes . Since the 
literature concerning this is very insufficient and many important as- 
pects either are not considered at all or are reported incorrectly, it 
seems that a short discussion is not superfluous. Especial consideration 
will be given to the various assumptions and neglections that are neces- 
sary in deriving the fundamental equation of Maxwell. The experimental 
work available, which is very insufficient and in part poorly dependable, 
can be used as an accurate check on the theory only in very few cases . 


I. THE FUNDAMENTAL EQUATION 


The theory of the evaporation process in a gas atmosphere owes its 
beginnings to Maxwell (reference l) and Stef sin (reference 2) . The theory 
rests on the assumption that the vapor in the immediate neighborhood of 
the liquid surface is completely saturated, and that consequently the 
velocity of evaporation depends simply on the velocity of. diffusion of the 
vapor into the surrounding space. This viewpoint was experimentally 

*The term ’Velocity" as used, in. this report corresponds to "rate” 
as currently used in American reports. 

1 From Physikalische Zeitschrift .der. Sowjetunion, vol. 6, 1934*. 
pp. 224-243. 
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substantiated by Stefan himself,, and also by Winkelmann’ (reference 3) and 
except for a limitation later discussed further, is regarded as an estab- 
lished fact. 


Since diffusion and heat conduction are physically related concepts 
and are also completely equivalent processes from a mathematical view- 
point, use can be made of the heat transfer; theory, which has been worked 
out in great detail, in the theoretical calculation of the velocity of 
evaporation. 


■.Now, turn to the derivation of the fundamental equation for the ve- 
locity of. evaporation of small droplets in a gas atmosphere. To do so, 
a series ..pfsimplifying assumptions must 'be made. - Later, these assump- 
tions will be examined, one * by one and corresponding corrections will be 
introduced in the fundamental equation. The assumptions are: 

1. The drop is" spherical. r 

.... 2. There is no motion of the drop relative to the gas atmosphere. 

3. The" atmosphere extends unbounded in all directions. 

4. The atmosphere is all at the same temperature and pressure. The 

lowering of temperature of the drop is neglected because of : 
the evaporation. 

5,. The evaporation process is stationary. 

6. The vapor is saturated on the surface of the drop. 


7 . The vapor pressure of the drop is vanishingly small in comparison 
with the total pressure. 


Since the process takes place sphero-symmetrically, Fick's law. 


|2 - Dio 
dt • 


expressed in spherical coordinates, becomes 

S(cr) _ ^ S £ (cr) 
St dr 2 

In the stationary instance this reduces to 


S 2 (cr) 

Sr 2 


= 0 


( 1 ) 


( 2 ) 


(3) 



NACA TM No. 1160 


3 


or, integrated 

0 * A + . - {k) 

If c G is the concentration of the saturated vapor, and Cx the concen- 
tration at infinite distance from the drop, then there is obtained: 


c - Cx * (c 0 - Ci) & 

r 


(5) 


where a is drop radius. 

The amount of substance (expressed in moles) diffusing away per unit 
time is (reference 1 ), 


« dc 

I Q = -knDr^ ~ <* 4naD(c 0 — Cx) 


(6) 


Therefore the velocity of evaporation of the drop in a gas atmosphere, 
as opposed to evaporation in a vacuum, is proportional to the diameter 
and not the surface. (Results are substantiated by Sresnewski (reference 
4) and Morse (reference 5*) Topley and Whytlaw-Gray (reference 6 ) 
checked experimentally formula ( 6 ) . Since 

r _ 7 dv 


where 

7 liquid density 

v drop volume 

M . molecular weight of liquid 
Then, from equations ( 6 ) and ( 7 ) 

aisfl . _ Si P S L . k , a constant, ( 8 ) 

dt 7 

that is, the drop surface changes linearly with time. Here, equation ( 6 ) 
is assumed to hold even for constantly changing drop Bize. This is true, 
down to a certain drop size, for the most widely different experimental 
conditions . 


„ d(a g ) 

M dt’ 


(7) 
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A correction must be applied to equation (6) if the previous assump— 


tions are not fulfilled. < Also, without loss of generality, 
written for ci. Then equation ( 5 ) becomes 

0 can be 

c - c 0 a/r 

(5’) 

and equation (6) 


I 0 = 4jtaDc 0 

(6«) 


II. FINITE VAPOR PRESSURE 


Firse, let fall the assumption that the vapor pressure of the liquid 
is vanishingly small compared to the gas pressure. Since the total pres- 
sure is everywhere the same, so must there he a concentration lowering of 
the gas equal and opposite to that of the vapor. Since, however, in the 
stationary state no flow of the gas can take place, so must the diffusion 
of the vapor he compensated by a convection current of the gas mixture 
directed toward the outside. The velocity of this current will be deter- 
mined by the following equation; 


C'U « D 1 


dc’ 

dr 


where the concentration of the gas has been designated as c*, and the 
diffusion coefficient of this in the vapor as D’ . But D’ = D and 


^ = _ do. Therefore 
dr dr 


u 


D dc 
c" dr 


(T f ) 


The total vapor flow then becomes 

I = -4nr 2 (d - cu) = -4«r a D f~ (^1 + (S') 

If C is designated as c + c 1 , the constant total molar concentration, 
then 

T _■ . 4nr a DC dc 
C - c dr 
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from which (reference 2) 

I « -hjtaDC In ( 1‘ ~ (9) 

(Concentration of saturated vapor is Co.) (For r = », c* = 0.) 
Expansion of the logarithm leads to 

I ? W*o (1 * - I„ (l * (10) 

As a first approximation, the percent errors involved In using (6* ) 
instead of (10)' is c o /'20. For water vapor at 20° C and atmospheric 

pressure, this amounts to only 1,2 percent. 


III. LIMITED SPACE OF EVAPORATION 


Here will he considered only the case of absorbing walls, on which 
there obtains a constant vapor pressure. With nonabsorbing walls the 
process is obviously not stationary and will be dealt with in another 
section.. . ....... 

If, for the sake of simplicity, it is assumed that the vessel is 
•‘Spherical and' that the drop is" located in the center, then there is ob- 
tained, Instead of (5) and (6), 


c = cp 


(c 0 - cp)a R 


r 


R — a 


( 11 ) 



4k(c 0 - cp)aD 




( 12 ) 


where R is the vessel— radius, and the constant vapor concentration at 
the wall is designated by cp. A vessel large with respect to the drop 

size has practically.no influence on the velocity of evaporation, 
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IV. THE LOWERING OF TEMPERATURE OF THE DROP 
BECAUSE OF ITS EVAPORATION 

This correction is the most important for volatile liquids. For the 
derivation of this, first of all, the radiation of heat and convection 
currents will "be neglected, and only the flow of heat to the drop "by con- 
duction will be considered. If it is also postulated that the coefficient 
of heat conduction is constant at all points of the gas atmosphere — that 
is, independent of the temperature or composition of the gas— vapor mixture 
within certain limits - then, in the stationary case the distribution of 
the temperature and the magnitude of the heat flow are expressed, respec- 
tively, in' the same formulas as the distribution of concentration and the 
velocity of evaporation, and can be expressed as: 

Ti - T.» &(?! - To) (13) 

Q = 4nKa (Ti - T 0 ) (l4) 


where T and Ti are the temperatures at r and at infinite distance, 
and T 0 is the drop temperature with K the coefficient of heat conduce 
tivity and Q the heat flow per unit time. The heat used in evaporation, 
however, is: 

Q« = I 0 1 = 4jtalDc 0 ( 15 ) 

(1 * molar heat of evaporation)' where D is assumed independent of tem- 
perature . 

(Translator’s Note: By virtue of its. evaporative potential, a drop 

initially at the same temperature as its surroundings, if quite volatile, 
will evaporatively self-cool to a temperature below that of its surround- 
ings such that the heat flow from the outside will supply the latent heat 
of vaporization. ) 

In the stationary state Q must equal Q,’ or (reference 1), 
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P 0 = f(T 0 ) (17) 


By (16) anl (17) T 0 and p D are clearly determined and are inde- 
pendent of the size of the drop. The above-mentioned linear relationship 
(equation (8) between the surface of the drop and the time remains undis- 
turbed by the lowering of the temperature. 

In order to coup are these formulas with experimental results, or to 
determine the value of the diffusion coefficient, (16) and (17) are not 
solved, but there is obtained directly from (Ik) , (15), and (7) 

Ti - To = (18) 

WKa 2MK dt 

there can be calculated Tj. — To, 
D (from (8)). 

For small T* — T 0 equations (l6) and (17) can be solved by means of 
the Claus ius-Clapeyron equation, and the following explicit expression for 
the temperature lowering can be obtained: 


From the measured values of 

dt 

T 0 and p Q (from (l'f) ), and, finally. 



(19) 


This correction was used by Topley and Whytlaw-Gray (reference 6) 

In their investigation of the rate of evaporation of small spheres of 
iodine by means of a spiral-spring balance.' The values of D' obtained 
in the manner described at various temperatures (14° to 30° C) differ 
only about 1 to 3 percent from values determined directly. The tempera- 
ture lowering of the small spheres in this investigation was, however,.; 
comparatively small (0.3° to O.5 0 C). 

Similar measurements were recently made by Houghton (reference 7) 
on water, droplets (a = — 300(i), which were hung on . very fine glass 

fibers. The decrease in volume of the drop was measured microscopically. 
The conventional psychrometric formula which Houghton used for calculat- 
ing the temperature lowering led, of course, to markedly varying values 
of D. Since, in this case, the temperature lowering has a significant 
value (about I5 0 at 25° C and in dry air), a ’study must be made "of the 
dependence of D on temperature.. As a first approximation, the as sump—., 
tion can be made that the coefficient of heat conductivity is constant, 
as is the distribution of temperature (13)/ and then the average value 
of 1/2 (Hi + K 0 ) can be taken for K in (16). The integration of equa— . 
tion (6) yields — for a variable D: 
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I • -Ks- (20) 

By differentiating (13)# there is obtained: 

dr = dT 
r 2 a(T x - To) 

For D it can be said, approximately. 



(Trans la. tor* a Note: Other investigators give 

D 0 Vt-/ . . 

See, for example, Sherwood, "Absorption and Extraction," McGraw-Hill 
Book Co., Inc., 1937*) 

By substitution of this expression in (20) there is obtained: 

I - (*) 

' L D 

(where Dj, corresponds to temperature T x ) . 

Therefore D must be replaced in equation (6) by the geometric mean 
/SjDo; but for small temperature lowering the ordinary arithmetic mean 
can be used. 

Use of that curve., will now be made in Houghton’s results, which 
shows the . evaporation of a water droplet at 21.7° C in perfectly dry air, 

and corresponds to a ■■■•- * 3-38 X 10 square centimeters per second. 

dt 

If-.^here is set, in equation (18) , 7 * 1, 1 * 10600 calories, 2M * 36 , 

K 13 = 6.0 x 10 cal/cm sec° . (ref erence 8) then. 
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T x - T 0 = I6.56 0 ; T 0 * 5.14° C ; c 0 = 3*79 x 10 7 moles per cubic centimeter 
from which, according to (8), where cj. = 0, 

D ■ A>oD ± = 0.248 

Finally, there can be calculated for . D at 0° C, according to (21), 
the value D =* 0.224, which is, a good approximation of the tabulated 
value 0.220 (reference 9). On the other hand, the curves given by Houghton 
for evaporation in partially Saturated air show strongly varying values of 
D; there could be considerable measurement error in the determination of 
the vapor concentration. 

In order to evaluate the magnitude of the heat transfer by radiation, 
which has been neglected, it will be assumed that the drop and also the 
walls of the vessel are perfect black bodies. Then the heat influx to the 
drop by radiation can be expressed by the Stefan— Boltzmann equation: 

Qx = 4jta 3 or (Tx 4 - T 0 4 ) 

If this equation is divided by (l4), then: 

cr ka ° T3 - 3 

Q K 

Ti 

(Translator's Note: This approximation assumes that T c < -jj~.) 

By substituting the values o = 5.7 x 10“ 5 and K = 2540 absolute units, 
(air at T x * 290° K), there is obtained: 

Qi 

-q = 2 r 19» 

Therefore, the neglecting of radiated heat in the calculation of the 
temperature lowering can introduce an error of only 2 percent at most for 
a = 100 4. According to an investigation of Langmuir (reference 10) the 
effect of thermal convection can also be neglected. 
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V. CONCENTRATION JUMP AT THE CROP SURFACE 


The unanswered question of Stefan and Winkelman concerning the exact 
value of the concentration of vapor at the surface of an evaporating body 
can be solved in the usual manner of gas kinetics. Quite analogous to 
the changes of velocity and temperature, respectively, present at a solid 
wall, there is a Jump in concentration upon the start of evaporation — 
which was first reported by Langmuir (reference 11) . This Jump at atmos- 
pheric pressure generally has a vanishing value on a flat surface, but 
it is quite otherwise in the case of small drops, when the diameters be- 
come comparable with the mean free path A of the vapor molecules. 

One of the molecules escaping from the surface will be, after travel- 
ing the distance A, at an average distance A from the surface. The 
ratio p = A/A certainly depends on the value of A/a and lies within 
the limits 2/3 (for a » A) and 1 (for a « A). Per unit time 
4ita 2 vc 0 molecules- escape from the drop surface, if vco is the number 
of molecules impinging on a square— centimeter-pen-unit time at a gas 
concentration Co • Here the usual presumption is made that the coeffi- 
cient of accommodation of the vapor molecules on the liquid surface is 
exactly 1. At the same time 4jta2vci vapor molecules condense on the 
same surface, where c x is the vapor concentration at a distance A 
from the surface. Therefore, the resultant amount evaporating is equal 
to Ii = 4na 2 vm(co — c*). This must equal the amount carried through by 
diffusion, I 2 = 4n(a + A)Dcj. . If there is substituted for D the ex- 
pression of Meyer 4/3 Avm, then there is obtained: 



•1 + 4/3A(a + pA) 
a 2 


(23) 


I * 


Io 


1 + 



PA 

a + PA 


(24) 


For the larger small drops 


(a .» A, p = | 


therefore: 


I = 


Io 


1 + 2* 
3 a 


(25) 
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.However, for the very email drops (a « A, 0*1) 

4jraDc 0 4jra 2 c 0 vm ^ 4jta 2 c 0 vm 


I » 


i ♦ s* - 


3 a a + A 


X + 2. — — ■■ a 

4 A 4 a + A 


1 + 1-sSL 


( 26 ) 


Therefore, very small drops evaporate in a gas atmosphere with the 
same velocity as in a vacuum. Moreover, in solving equation (24) it is 
possible also to start from the velocity of evaporation in a vacuum, in 
that this is multiplied, by the probability that an evaporated, molecule 
does not again encounter the drop in its subsequent motion. 

Equation (25) shows that the usual linear change of drop surface 
holding for small drops is no longer valid. Since A at atmospheric 
pressure is of the order of 10" 5 centimeter, the correction calculated 
by (25) is already 7 percent for a = lp, and 40 percent for a = O.lp. 

A decrease in the rate of change of the drop surface, which began to be 
noticeable at a = lp, was actually observed by Speakman and Sever 
(reference 12) on a series of organic compounds . The explanation by 
these authors that this phenomenon was due to a lowering of the vapor 
pressure of the drop by dissolved nonvolatile impurities, can perhaps 
also be true. Unfortunately, the curves presented in the above-mentioned 
book do not permit a quantitative comparison with equation (25). 

In a recent work by Woodland and Mack (reference 13) an opposite 
observation was made — an increase in d(a 2 )/dt upon diminishing the 
size of the drop. The cause of this discrepancy must lie in something 
else; in the proposed explanation of their findings, however, these 
authors come to very erroneous conclusions, such as the presence of a 
layer of vapor, 0.5p thick, surrounding the droplet, , and so on. 

It is easy to see that the correction of the change of vapor pres- 
sure as a consequence of the curvature of the., surface and the charge of 
the drop can be neglected in comparison with the correction Just discussed. 


VI . NONSTATIONAHI PROCESS 


The exact calculation of the velocity of evaporation in the noneta— 
tionary case meets, in general, the; greatest mathematical difficulties. 
Fortunately, in practice most of the conditions of the process can be re- 
garded as quasi-s tat ionary (that is, the process has at every instant 
the same velocity as in the stationary state which corresponds to the 
boundary conditions at that instant) as follows: 
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Into the infinitely extended gas atmosphere, in which .the vapor con- 
centration is everywhere zero, a drop of radius a is introduced at time 
t = 0. Apart from the .Jump in concentration; at the surface, and so forth, 
the primary problem is the solution of the differential equation (2) with 
the following boundary and initial conditions, respectively: 

c = 0 at t = 0 and r > a, 
c * c 0 at t > 0 and r a a. 


The solution is given by the following expression: 


2c pa 
r Jn 


e dx 


r - a 

2 /St‘ 


From' this there is obtained. 


I = - kite . 2 D — = httaDcp (\ + •—=) 
dr \ /nSt/ 


dr 

(r=a) 


» I 0 (l + = 

\ /fiDt, 


Since ^nL at atmospheric pressure is of the order of magnitude 1, 1 

~ f a\ 

there can be substituted I = I n i 1 + — /. The correction amounts to 

only 1 percent after 1 second, even in a heavy fog (a = IOO4) . In order 
to decide to what extent the evaporation of a drop: can be regarded as 
stationary, the' time interval t x is compared, • after which the correction 

term a//«Dt reaches a definite small value A, with the time t 2 
necessary for complete evaporation of the drop. 


tv. - 


A 2 nD 


and from (8) 


2MDc 0 
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therefore 


• ■ t 2 'A«Dr ' ' . ‘ . 

If there are substituted the values A * 0.01/ M » 18, and 1 

_7 • 

c D *= 3-79 x 10 moles per cubic centimeter, vhich correspond to the 
evaporation of vater drops in dry air at 21.7 C (see above), then 

= 0.043, or in other words, even after the course of the twentieth 

part of the time of total evaporation the velocity of evaporation exceeds . 
the velocity corresponding to the stationary state by only 1 percent. In 
damp air or with less volatile liquids the approximation to the stationary 
state is realized still more quickly. The effect of the gradual reduction 
of drop size which was left unconsidered in the discussion will be de- 
scribed further. 

For the calculation of the velocity of evaporation of droplets of 
solutions, the nature of the variation of vapor pressure at the surface 
with time beoomes important. If it is assumed that the concentration of 
the saturated vapor is a time function c(t), then, instead of (27), 
there is obtained 


t — e x dx (reference 15) 

ilThr 2 - 


r — a 

2 ✓Tit 


which can be easily verified. From this is obtained (omitting the repe- 
tition of the somewhat detailed calculations): 


c(0) + at 


I * -4na s D ~ - 4naDc(t) j 1 + 

or l- 

(r=a) 


^ e’[t(l - 


x s )]dx 


— (3D 


✓Dtit 


where ■■'C* means - y i 


*.4.1- '« V , »4 i'\V ■ M 

* .,.1 ‘ f "#ii; ‘ ’m ’ 5 "*■ 
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Whether or not an approximation to the stationary state is possible 
here depends on the special form of c(t). In general, the following can 
be saidj wherein it is necessary to take into account the most important 
practical cases of the decrease in c(t) with time. In order to rea ch a 
good approximation of the stationary state, it is sufficient that a//Dtn 

Increase by a definite amount A after a time interval t, which is 

small relative to the interval t 3 , meanwhile the factor times a/'/Stit 
remains of the order of magnitude of 1, If |c'(t)| is designated by 
f(t) and it is noted that: 


c(t) = c(0) 


+ i 


c* (t)dt = c(0) 


1 

— t / f (tx)dx 
J o 


then this factor takes the following form: 

i 


F = 


c(0) - 2t £ f [t(l - x^Jdx 

c(0) - t / f (tx)dx 

J o 


With decreasing f (t) 

.1 . l 


j' f(tx)dx = J' f[t(l-x)]dx> j' f[t(l - x 2 )]dx 


and for t 3 the interval will be taken after which c(t) is half of its 
original value, that is, the root of the equation 

= t f f (tx)dx 

In this interval F < 2. 

With increasing f (t), that value of t can be chosen for t 3 sudi 
that F vanishes; c(t) will then not reach l/2c(0). In this interval 
F < 1. It is sometimes more convenient to use t 3 ' < t 3 instead of t 3 ; 
so that the equation c(0) = 2tf(t) Ib satisfied. 

Now, consider the continuous diminishment of the evaporating droplet 
due to disturbance of stationariness . The strict treatment of this prob*- 
lem is extremely complicated; as an approximate estimate use can be made of 
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the following consideration. If there is originally a stationary state 
corresponding to drop radius a and vapor concentration c Q , then the 

states arising from a decrease in a of Aa at constant Co, and a de- 
crease in Cq of Ac 0 at constant a, respectively, are given by the 
following expression: 


* i {°° a - °<^ L 1 " * (■"'—)} 


and, respectively. 


c 



SAC 0 



= Gauss' probability function.) 

If, now, aAc 0 * CoAa, then these expressions can differ from each 
other only by an Infinitely small amount in the second order. If, there- 
fore, the time function c 0 (t) * c o (0)qp(t) is a good approximation of 
the stationary state, then the equation a(t) = a(0)cp(t) should apply 
for decrease in drop size with time. For the evaporation of the drop 
there can be substituted, as a first approximation, from (8) 

a(t) = </a(0) 2 — kt, Here then: 


f(t) 


da k 

dt 2/a(0) 2 - kt 


and for t 3 * there is obtained 


from which 


a(0) = 


t 3 'k 


A(0) 2 - kt*,' 


t » = 0.6 = 0.6 

3 -k 


^2 


(t 2 is time of total evaporation.) 

Then, return to a numerical factor 0.6 for equation (29); it follows 
from this, that there can be assumed a quasi— stationary course of evapo- 
ration of water droplets as a good approximation. This conclusion obvi- 
ously agrees with the experimental results of Houghton; - 
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Thus far, the evaporation of a single drop in an infinitely extended 
medium has "been considered. In practice it is necessary to do, however, 
with a large number of droplets. In this case the evaporation of a drop 
takes place as though it were in a vessel with nonabsorbing walls, the ' 
volume of which vessel being equal to the average volume enclosing each 
drop. The rigorous solution of this problem is very difficult. Practi- 
cally, for the small weight concentrations occurring in fog, that is, a 
large mean distance between the drops relative to their diameter, a sim- 
ple approximation will serve. Then, assume that the droplet is in the 
midpoint of a spherical vessel of radius B with nonabsorbing walls. Tte 
distribution of concentration for different times is shown in figure 1 
with dashed curves; the solid curves correspond to the stationary state 
(for absorbing walls with various vapor pressure, see sec. III). 
(Translator's Note: The solid curve represents the stationary state for 

the same boundary conditions that exist at a particular instant in the 
nonstationary state.) It Is easy to see that a curve of the first type 
must. In fact, lie under the curve of the second type with the same end- 
points. There results, therefore, for the velocity of evaporation in a 
known time t + At, a lower limit. (See equation (12).) 


4«aD(c 0 - c t +&t) 
lx = 


An upper limit can be obtained as follows. If the drop is placed 
at time t in an unbounded space with vapor concentration ct, then its 
evaporation will take place more quickly than in the vessel. In time 
t + At it will reach the following value: 


Is 


= 4naD(co — ct)l 


1 + 




(see equation (28)). 

If the time interval At is chosen- such that 


i - - 

B 

then the pictured curves are actually obtained, which do not intersect 
each other more than once. The equivalence of the two expressions for 
the amount of material evaporated in time interval At and stored in 
the gas atmosphere, respectively, yields: 
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1? 


B 

^naL(c 0 - c t )At ~ kit (ct + ^t( r ) ~ c t( y ) 
1 - - • : a k 

B 


r e dr 


where ct(r) is the value ct corresponding to the stationary distribu- 
tion. By substituting the expression for c(r) (equation (11)), there 
is obtained 


Ac * C t+At c t 


3Da(c 0 - c t )At 

— V- 5 — 


(32) 


or, since 


At 


Bf 

nD 


Ac 


(°0 C t^ R 


(33) 


At small values of a/E the stationary curves give a good approx- 
imation. With the same approximation (32) can be regarded as a differ- 
ential equation. By integration, obtain 

-3Dat/B 3 -3Dat/B 3 ' ' 

I = M-jtaDco© - Io© 


This equation is obviously only applicable when the decrease in vol- 
ume of the drop can be neglected (at very small vapor pressures or pres- 
sure differences). Otherwise, the system of differential equations ..(8) 
and (32) must be solved, which usually presents no difficulties. 


VII. MOTION OF TBE DROPLETS 


The exact calculation of the velocity of evaporation of a droplet 
moving with respect to the gas medium is scarcely possible. In a coordi- 
nate system relative to the drop (that is, drop at r * 0, Z = 0, a = 0), 
this problem leads to the differential equation: 

” = I&c — v grad c (35) 

at 

where with v is meant the velocity vector of the gas stream, the value 
of this vector being determined in space by the nature of the flow. 
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Even in the simplest case of laminar flow, for which Stokes 1 law holds, 
and a stationary state, equation (35) is practically insolvahle. It can 
he solved only for an infinitely small flow velocity: if this velocity 
is called V at a great distance from the drop, then let 


c 


c o a 

r 


+ V<p 


where the term Yep represents the disturbance of the concentration dis- 
tribution caused by the flow. From (35) there is obtained (with 
bc/bt = 0): 

DV£cp + v ( + V grad cp } = 0 

\ r / 

and, since the last term can be neglected; 


- 


CqQ, — 


Dr 3 Y 


vr 


If the Stokes expression for v is substituted in- this equation, 
the following result is obtained: At two points which are symmetrical 

with respect to the plane passing through the drop middle perpendicular 
to the flow direction, £cp possesses the same absolute value, but re- 
versed sign. Since further, cp vanishes on the surface of the drop and 
at infinite, distance from it, then the application of the potential theory 
leads to the conclusion that at two diametrically opposite points on the 
surface grad 9 has the same length and direction, that is, that the 
acceleration of the evaporation produced by the stream on one side of the 
drop is exactly compensated by the slowing down of the process , on the 
other side. 


This conclusion could approximately hold for a flow velocity, wherein 
the transfer of material by convection can be neglected in comparison 
with the transfer by diffusion — at least at distances not too far from 
the drop, that is, when cjv| « d| grad c|. If by substitution here 
c = coa/r, there is obtained. 


c 0 aV 

r 


« 


Dc 0 a 

1.2 ’ 


aV 

or — « 1 
D 


(a - r) 
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Since it is known that D = ^ (t) is coefficient viscosity of the 

medium, p its density), then aV/D becomes the familiar Reynolds 
expression apV/r\. For small Reynolds numbers (for example, for free 
fall of fog droplets) a vanishingly small influence of the motion on the 
evaporation velocity would be expected. The validity of this conclusion 
obviously can only be decided by an experimental procedure. 

In conclusion, it should be said again that all the considerations 
and conclusions carried out here are also directly applicable to the 
reverse process — the growth of drops in a supersaturated vapor. 


Translation by John Kelson Howard and Mitchell Gilbert, 
National Advisory Committee 
for Aeronautics. 
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